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Abstract 

We apply the momentum integrated Boltzmann transport equation to study 
the time evolution of various quark flavours in the central region of ultra-relativistic 
heavy ion collisions. The effects of thermal masses for quarks and gluons are incor- 
porated to take into account the in-medium properties of these ingredients of the 
putative quark gluon plasma. We find that even under very optimistic conditions, 
complete chemical equilibration in the quark gluon plasma appears unlikely. 
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Quantum Chromodynamics (QCD), the theory of strong interactions, predicts 
that at very high temperature and/or density, the colourless hadronic matter dis- 
solves into its coloured constituents, the quarks and gluons, such that the bulk 
properties of strongly interacting matter are governed by these degrees of freedom. 
Such a locally colour-deconfined phase of matter is known as Quark Gluon Plasma 
(QGP). It is expected that the temperature and density achievable in ultrarelativis- 
tic collisions of heavy ions is favourable for QGP formation, although transiently. 

Many aspects of this transition, e.g., the order of the phase transition, the value 
of the critical temperature, signals of the transition, thermodynamic equilibration 
, are still uncertain. So far, most of the works assume thermodynamic equilibrium 
throughout the evolution history, after an initial proper time Tj (~ 1 fm/c). Thus 
the thermodynamic quantities like pressure, entropy and temperature have direct 
meaning with respect to an equilibrated system. It has however been recently shown 
that the approach to equilibrium in ultrarelativistic heavy ion collisions proceeds 
through a succession of many time scales 1 ' 2 . In particular, the dominance of the gg 
cross section over the qg or qq cross sections 3 was argued to imply that the gluons 
equilibrate among themsleves appreciably earlier than the whole system of quarks 
and gluons. In a recent work 1 , we showed that the time scale of kinetic equilibra- 
tion for the light quarks (u, d) is ~ 1 fm/c; for massive quarks, it increases with 
the mass to rather large values. The significance of these considerations for QGP 
diagnostics has been discussed in references 4 and 5; the emission of particles from 
the pre-equilibrium era (r 9 < r < T t h), Tth denoting the time for full thermodynamic 
equilibrium, may indeed populate the invariant mass or p T windows thought to be 
relevant for signals from the thermalised QGP. 

In this letter, we study the chemical evolution of the quarks in the system 
from the epoch r g ( the proper time when gluons thermalise) onward. Chemical 
equilibration has very important implications for the signal of QGP formation 4,5 ' 6 
. It would be ideal to study this problem in the QCD based kinetic theory and 
attempts along these lines have been made recently 7 . Such calculations, appealing 
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as they are, rely heavily on rather involved numerical simulations and also suffer 
from dependence on parameters needed to mimic the non-perturbative effects. We 
propose, as in ref. 1 , a scenario which retains the essence of the kinetic theory to 
the following extent : the gluons thermalise completely at a proper time r g earlier 
than the quarks. The gluons carry about half of the total momentum and the quarks 
carry a tiny fraction. We restrict ourselves to the central rapidity region where the 
number of valence quarks is assumed to be negligibly small. We shall, however, 
relax this condition when we look at the situation at SPS energies. 

In order to treat the equilibration of the species one must follow the micro- 
scopic evolution of the phase space distribution function governed by the 
Boltzmann transport equation : 



noting a quark flavour. The Boltzmann equation is a non-linear, integro-differential 
equation for the phase space distribution function of the particles. For our present 
purpose, however, the problem can be addressed in a rather simple manner without 
losing much information on the time scales for chemical equilibrium. 

The decoupling of the relic particles in the early universe was studied 8,9 by 
integrating the Boltzmann equation over the momentum of the particle to obtain 
an equation for the evolution of number density. We use a similar approach to look 
at the number density evolution of the partons in ultrarelativistic collisions of heavy 
ions; the major difference between the two cases is that the expansion dynamics in 
the early universe (big bang) is governed by the Hubble expansion 9 whereas in the 
case of heavy ion collision (mini bang), the expansion dynamics is assumed to be 
governed by the boost invariant (Bjorken) solution of relativistic hydrodynamics 10 
. The equation for the number density evolution then reads 



where n q is non-equilibrium quark density, r the proper time and g q the statistical 
degeneracy of the quarks of flavour q. 



L{f q } = C{f q } 



(1) 



where L is the Liouville operator and C the collision operator, the subscript q de- 




(2) 
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Clearly, the time evolution of n q would be largely governed by the various pro- 
cesses determining the collision term in eq. (2). The most important contribution 
to this is from gluon fusion from {gg — > qq). 

Since we have considered complete equilibration for gluons, reactions like gg <-> 
ggg need not be included 1 ; the gluon number density at each instant is uniquely 
determined by the temperature of the system. However the gluons in a finite tem- 
perature system acquire a sizable thermal mass may qualify as real excitations 11 ' 12 . 
Thus the gluon decay process (g — > qq) and its inverse is expected to play a major 
role; in particular it has been argued in the literature 11 that for small quark masses 
(m e ff/T « 1), the gluon decay dominates over the fusion process. We must there- 
fore take the gluon decay into account while evaluating the collision term. We ignore 
the quantum effects as the temperature is expected to be very high but include rel- 
ativists effects for the same reason. Thus gluon are described by the relativistic 
Maxwell-Boltzman statisics, f g (E) = g g /(2ir) 3 e~ E l T . Then eq.(2) can be written 
as: 

eq 

where n eq (T) = g q /(2n) 3 J d 3 p/(e E / T + 1), is the equilibrium density. Eq.(3) is a 
particular form of the Riccati equation; there is no general closed form solution of 
this equation. 

In deriving eq.(3) we have assumed T (or CP) invariance, i.e., \M\ gg ^ qq = 
\M\ qq ^ gg and \M\ g ^ qq = \M\ qq ^ g . R gg ^ qq (T) is the quark production rate per 
unit four volume by the reaction gg — > qq and R g ^ qq gives the same for the gluon 
decay 11,12 , g — > qq . The thermal quark masses have been taken into account using 
m e q ff(T) = yj rricurrent + vp? th in all the rates. In the present context in a chemically 
non-equilibrated scenario m th is given by 13 

r o 2 T 2 

-4 = (i + ^)^n ^ 

where g is the colour charge, r q = n q /n eq . The full derivation of eq.(3) has been 
omitted here for the sake of brevity; it will be reported elsewhere 14 . Eq.(3) is 
a very convenient form of the Boltzmann transport equation with all the relevant 
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features for our present purpose. The first term on the right hand side of the 
equation represents the "dilution" of the density due to one dimensional expansion, 
the destruction of the quarks is proportional to the annihilation rates of qq and 
destruction is balanced by creation through the inverse reaction when n q = n eq . The 
last term represents the change in quark density due to gluon decay and its inverse 
process, arising from the thermal gluon mass. Note that there is no overcounting 
as the phase space integrations involved in the gluon fusion and the gluon decay 
channels are over different domains. For a detailed discussion see ref. 11. 

To solve eq.(3) for n q (r) we need to know T(t), which is obtained by solving 
the hydrodynamic equation 

£ + — = (5) 
ar r 

where P is the pressure and e the energy density. An equation of state relating P 
and e can be parametrized in the form 

P = cle (6) 

Then, for the non-equilibrium configuration, we can write 

e=[ar g + br q (T)]T± (7) 

where a = (8tt 2 )/A5cJ 2 and b = (7n 2 )/60 Nf c~ 2 . The medium effects enter eq. (7) 
through the thermal mass of the quarks; the higher the temperature, the higher the 
thermal mass. Therefore the equation of state may be parametrized as P = c 2 (T)e, 
where c; 2 = + 3] with g eff (T) = (e + P)/(2tt 2 T 4 )/45. This would 

amount to making a and b in eq. (7) functions of T. For the cases considered here, 
such effects are negligibly small. 

In chemical equilibrium, r q for all species should be unity and in correspondence 
with our premise, we take r g —1 all through. Nf is the effective number of massless 
flavours; appearance of quark masses amounts to reducing the actual number of 
flavours 15 . For the present, we confine our attention only to the three lightest 
flavours, u, d and s. The difference between the u(d) quarks and s quarks arises 
only in terms of their current mass. We have verified that because of the high initial 
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temperature and the corresponding large thermal masses, this amounts to a very 
small change in the effective Nf. 

Equations (5) and (7) together yield, 



T T 



a + br q (r g ) 



3/4 



T*r g (8) 



a + br q {r) 

which is the cooling law in the present non-equilibrium scenario. 

To solve eq.(3), we need to specify the initial conditions. There is no unique pre- 
scription for evaluating them at the present state of the art. Biro et al. estimated 16 
the initial parton density from the HIJING algorithm 17 . On the other hand, Shuryak 
and Xiong 5 , advocated, as in ref. 1, the use of "perturbative" QCD(pQCD), i.e., 
the deep inelastic structure functions for the estimation of the initial quark density. 
Given the importance of the issue, we have, in the present work, used both these 
estimates. 

In the pQCD approach, the initial values of n 9 (r 9 )'s are obtained from the 
integration of the structure function 18 , 

fx ■ uJx)dx 

n q = 2Ax JXmm qK — (9) 
7iR A r g 

where A is the mass number of the colliding nuclei (208 in our case for Pb), x m i n 
is the minimum value of x, taken 1 ' 5 to be 0.02, dictated by the applicability of 
perturbative QCD. It is well known that the structure functions depend on x and 
Q 2 and therefore, it is important to choose the proper range of Q 2 . The range 
of Q 2 appropriate for our present context can be determined by comparing the 
values of the running coupling constant a s (Q 2 ) with a s (T = T g ) 19 . Thus for initial 
temperatures T g of 500 and 660 MeV (corresponding to r g = 0.3 and 0.25 fm/c) at 
RHIC and LHC energies 2 , respectively, the appropriate ranges of Q 2 are 25 and 55 
GeV 2 . Then, the values of n q (r g ) estimated from eq. (9) come out to be 1.26 and 
1.33 fm~ 3 at RHIC and LHC energies. From ref. 16, these values are 0.7 and 2.8 
fm~ 3 , respectively. 

Eqs.(3), (4) and (8) were solved self consistently by the Runge -Kutta method to 
obtain the non-equilibrium density n q (r). As is evident from the above discussions, 
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the non-equilibrium density n q has an explicit dependence on r and an implicit 
dependence on r through T(t). But the equilibrium density n eg (T) has only an 
implicit dependence on r through T(r). The ratio r q thus assumes an universal 
feature, since the implicit time dependence gets eliminated. The time dependence 
of the ratio r q can then be used as a ready marker for chemical equilibrium; the 
time at which the explicit time dependence of r q vanishes, simultaneously with r q — > 
1, corresponds to the time for chemical equilibration for the flavour q. 

In contrast to the earlier work of Biro et al 16 , we have not included reactions 
like qq — > QQ, since the initial system is dominated by gluons 2 ' 7 ; the quark density 
is very low compared to the gluons. These reactions are suppressed by a factor of 
about 1/16 compared to the case when the quarks are in complete equilibrium 5 . 
Furthermore, we would like to emphasize at this juncture that we are looking at an 
optimistic scenario, where the gluons have completely equilibrated already. 

The relative importance of the gluon fusion (gg — ► qq) and gluon decay (g — > qq) 
channels is obviously a crucial issue in the present context. In fig. 1, we show 
the reaction rates for these two channels. If one considers the thermal contribu- 
tion to the quark masses also, then the effective quark mass becomes m e q ^ ' (T) = 
\Jm1 urrent + mf h . One can readily see from fig. 1 that the gluon fusion dominates 
over the gluon decay modes for the entire range of interest; it has however to be 
noticed that the higher the m e '* the greater the difference between them, implying 
that for light quarks(u,d and even s), the decay channel does play an important role. 

The ratios r q = n q /n eq are plotted in fig. 2 (a) for RHIC energies for the initial 
conditions specified by eq. (9). At these energies, the ratios for u(d) quarks and s 
quarks are the same at early times but at late times the u(d) quarks dominate over 
the s quarks. The reason is that at early times the thermal mass (oc T) dominates 
over the current quark mass for all flavours. At later times, however, the thermal 
mass decreases with decreasing temperature and the difference in the current masses 
of u(d) (~ lOMeV) and s (~ lbOMeV) quarks becomes important. There occurs 
thus some Boltzmann suppression for s quarks relative to u(d) quarks. 
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The same calculation has also been carried out for the initial condition specified 
by the algorithm HIJING. At RHIC, the initial r q for all three flavours are lower 
than those estimated from eq.(9) (fig. 2(a)). As is to be expected, this makes 
equilibrium even less likely, which is borne out in our detailed calculation. For 
LHC, the initial r g 's are systematically higher than those obtained from eq.(9); one 
could thus naively expect that equilibration may become easier. We therefore show 
in fig. 2 (b) the result for this case also. One sees no sign of saturation with r for 
either u(d) or s quarks. 

In fact, one can also estimate the initial quark density from the prescription 
N q = o q ■ ^4^(0) where 7^(0) is the overlap function for central AA collision and 
<j q the production cross section for the quark 20 q , one would obtain significantly 
lower values of n q (r g ), or equivalently r q (r g ). {e.g. n q (T g ) ~ 0.48 at RHIC compared 
to 1.26 obtained from eq.(9)). This would then imply lesser probability for chemical 
equilibrium. We have verified that this is indeed so; the qualitative behaviour of 
the results shown above remain unaltered while quantitatively the conclusions are 
strengthened at all the energy domains considered here. 

We have so far considered futuristic scenarios at RHIC and LHC where the 
assumption of a baryon free region is perhaps applicable. For the presnt, namely 
the SPS energies, such an approximation is not at all appropriate; it is expected 
that there would be a substantial amount of stopping at these energies 21 . Never- 
theless, study of flavour equilibration at SPS energies becomes a very pertinent and 
timely issue. We have therefore investigated two extreme scenarios- one of total 
transparency as above and one of complete stopping. 

In fig. 3 we show the time evolution of r q and T as a function of r for the 
transparent case. The initial conditions are determined by eq.(9) alone, as, to our 
knowledge, the corresponding estimates from HIJING do not exist in the literature. 
We find that at early times r s is slightly larger than r u or r d ; this is due to the initial 
normalisation. The non-equilibrium density n u ^) ~ n s but the equilibrium density 
for s quarks is less than that of u or d quarks because of the difference in their 
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masses. This results in r s being larger than r u ^)- Expectedly, this effect washes 
out with progress of time. In any case, there is no hint of flavour equilibration at 
SPS energies either, under the assumption of complete transparency. In the other 
extreme case, that of complete stopping , however, the situation is most interesting. 
Because of the presence of a large number of valence u and d quarks, the initial 
values of r u and r d are very large. For central collisions of lead nuclei, the initial 
value of r„(d) turns out to be 0.98(1.06), taking into account both valence and sea 
contributions. As r increases, r q tends to 1, r u from below and from above. 
Nevertheless, the stationarity of r q with r does not seem to be achievable within 
the life time of the QGP. It is thus fair to conclude that at SPS, u and d quarks, 
although not in equilibrium, may not be too far from equilibrium in the event of 
complete stopping. The situation for the s quark is, however, quite different as all 
the initial s quarks come from the sea. Thus r s (r g ) is still rather small. But now, 
one has to include uu(dd) — > ss, in addition to gluon fusion and decay channels. 
We observe that even for such a favourable situation the time evolution of r s is 
indistinguishable from that shown in fig. 3. 

These findings are obviously in contrast with the earlier expectations of other 
authors 22,23 where s quarks were expected to come to equilibrium at SPS energies. 
The difference is ascribable to the effect of the thermal masses, which must be 
included if proper account of the in-medium effects is to be taken. We have shown 
here that these effects do play a very important role indeed and may affect the 
conclusions nontrivially. Obviously, implications of these issues in QGP diagnostics 
is an urgent task; work along these lines is in progress. 
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Figure Captions 

Fig. 1 : Comparison of quark production rates due to gluon fusion and gluon decay, 
the value of a s is calculated from the parametrisation of ref 19. 

Fig. 2 : Ratio of non-equilibrium density to equilibrium density and temperature as 
functions of time (a) at RHIC energy for initial condition given by eq.(9) and (b) 
at LHC energy for intial condition given by eq.(9) and ref. 16. 

Fig. 3 : Same as fig. 2 (a), at SPS energy (see text). 
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